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A TIME-DOMAIN METHOD FOR DETERMINING THE DEPTH 
OF A FLUOROPHORE IN A TURBID MEDIUM AND 
ESTIMATING THE FLUROPHORE CONCENTRATION 

CROSS -REFERENCE TO RELATED APPLICATIONS 

[0001] This is the first application filed for the present 
invention. 

TECHNICAL FIELD 

[0002] . This application relates to fluorescence 
measurements in turbid media and more specifically to 
fluorescence measurements using the time domain optical 
method. 

BACKGROUND OF THE INVENTION 

[0003] Optical fluorescence imaging of turbid media, e.g. 
biological tissue, is primarily achieved with Continuous 
Wave (CW) methods. Typically a light source is employed to 
illuminate the object of interest, e.g. a mouse in vivo, 
and the emitted CW fluorescence intensity signal from the 
flurophore is measured directly with a CCD camera, e.g. 
AntiCancer Inc. 

[0004] In this context, it is often useful to determine 
the concentration and depth of the fluorophore. However, to 
assume that the direct CW fluorescence intensity signal is 
proportional to the flurophore concentration can be 
misleading since the depth of the flurophore will also 
impact the. CW fluorescence intensity signal. In fact, given 
a single CW source and CW detector measurement it is 
impossible to decouple flurophore concentration and depth. 



BEST AVAILABLE 




10/31/2003 16:09 FAX 613 230 6706 



OG1LVY RENAULT 



@014 



File No. 15186-46USPR 

- 2 - 



[0005] To address this problem of decoupling depth and 
concentration in CW measurements, some researchers are 
exploring CW tomography, in order to reconstruct a three- 
dimensional distribution of the flurophore concentration in 
the object, thereby attempting to account for the depth of 
the flurophore. However, these tomographic approaches 
require multiple source-detector pair measurements from 
many angles, combined with complex, computer intensive 
inversion algorithms. Furthermore, CW tomography requires 
an assumption about the scattering coefficient of the 
object, since CW can neither decouple the intrinsic 
absorption and scattering of the object. 

[0006] The most complete description of photon migration 
in turbid media is provided by Time Domain (TD) optical 
methods which have previously been used to decouple the 
attenuation coefficient, given from CW intensity 
measurements, into the underlying absorption and scattering 
coefficients. However TD methods have not been applied to 
obtain depth and concentration of fluorophores. 



SUMMARY OF THE INVENTION 

[0007] The invention relates to the measurement of 
fluorophores position and concentration within a turbid 
medium. More specifically the invention relates to the 
measurement of fluorophores position and concentration 
within a turbid medium using Time Domain (TD) optical 
method. 

[0008] There is provided a method wherein TD optical 
methods are used to decouple fluorophore depth and 
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concentration from a single source and detector measurement 
employing a direct analysis technique, thereby avoiding the 
aforementioned burden of CW tomography and overcoming the 
limitations of conventional CW direct fluorescence 
intensity imaging. 

[0009] In one embodiment there is provided method for 
determining the depth of a volume comprising a fluorophore 
in a turbid medium using time domain (TD) optical 
fluorescence by obtaining Temporal Point Spread Function 
(TPSF) , data corresponding to at least one TPSF, by 
injecting light at an injection point at an excitation 
wavelength of the fluorophore and detecting light at a 
detection point at an emission wavelength of the 
fluorophore, the injection and detection points being in a 
reflection geometry and substantially equidistant from said 
volume, determining a time (t ma2C ) substantially 
corresponding to the maximum of the TPSF and correlating 
tmax with the depth. 

[0010] In another embodiment there is also provided a 
method for estimating the concentration of a fluorophore in 
a volume of a turbid medium using optical fluorescence by 
obtaining depth of the volume, providing optical properties 
for the medium, obtaining a CW intensity surface reflection 
measurement of the fluorophore, and normalizing the CW 
intensity measurements using optical properties of the 
medium and the fluorophore depth to obtain a relative 
fluorophore concentration and/ or absolute concentration. 

[0011] in yet a further embodiment there is provided a 
method for generating a tomographic image of a fluorophore 
distribution in a turbid medium by obtaining a topographic 
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image of the fluorophore distribution, determining the 
depth of a plurality of volumes of interest comprising the 
fluorophore using depch determination methods of the 
present invention, and combining the depth information and 
the topographic image to generate a tomographic image of 
the fluorophore distribution. 

BRIEF DESCRIPTION OP THE DRAWINGS 

[0012] Further features and advantages of the present 
invention will become apparent from the following detailed 
description, taken in combination with the appended 
drawings, in which: 

[0013] Pig- 1 is a schematic representation of a preferred 
fluorescence signal acquisition; 

[0014] Fig. 2 is a 3 -dimensional graphic relating t maJC / 
depth of fluorophore and the concentration of the 
fluorophore ; and 

[0015] Fig. 3 is a 3 -dimensional graphic relating CW 
signal intensity, depth of fluorophore and the 
concentration of the fluorophore. 

[0016] It will be noted that throughout the appended 
drawings, like features are identified by like reference 
numerals . 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
[0017] In the following description it will be appreciated 
that CW intensity can be obtained using direct CW 
measurements or by integrating a Temporal Point Spread 
Function (TPSF) obtained by the time domain (TD) method. 
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[0018] It will also be appreciated that the term 
fluorophore can refer to either an extrinsic fluorophore 
which is understood to be a fluorophore Chat is added to a 
given medium or an intrinsic fluorophore which is 
understood to be a fluorophore that is normally comprised 
in a given medium. For example biological tissue may 
comprise molecules that naturally fluoresce and are 
therefore intrinsic f luorophoreB . 

[0019] The fluorescence emission from a fluorophore 
embedded in a turbid medium can be measured over time using 
time domain optical fluorescence to generate Temporal Point 
Spread Function (TPSF) data. Figure 1 depicts a preferred 
arrangement used to obtained the TPSF. One or more 
fluorophore is included in a turbid medium and excitation 
light is injected at injection point 10. The photons 
diffuse and a certain fraction of the photons eventually 
reach the fluorophore molecules contained in a volume of 
interest (VOI) 14 which are thereby excited. The molecules 
then emit a fluorescence signal. The emission photons will 
diffuse and a fraction of them will reach detection point 
12 to produce TPSF data. 

[0020] TD fluorescence experiments were conducted with a 
single source -detector measurement for a small inclusion of 
fixed flurophore concentration submerged into a turbid, 
liquid medium. It was observed that the temporal position 
of the maximum of the TPSF generated by the fluorescent 
photons (tnax) increased with the depth of the submersion of 
the inclusion in the turbid medium (see figure 1) . 

[0021] TD fluorescence experiments were also conducted 
with a single source- detector measurement for small 
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inclusions of various fluorophore concentrations submerged 
at fixed depths in a turbid, liquid medium. It was observed 
that t,,^ generated by the fluorescent photons remained 
constant for inclusions of the same submersion depth, 
regardless of their fluorophore concentration (also shown 
in Figure 2) . Hence it was demonstrated that t^ from a 
fluorophore inclusion permits decoupling of its depth from 
its concentration. 

[0022] For the above series of experiments CW data was 
also generated by temporally integrating the complete 
fluorescent TPSF. As expected, the CW fluorescence 
intensity decreased as the depth of the fixed fluorophore 
concentration increased, but it also decreased as the 
fluorophore concentration of inclusions at a fixed depth 
decreased, as shown in Figure 3. Hence, it was demonstrated 
that direct CW fluorescence intensity measurements alone 
cannot decouple fluorophore inclusion concentration from 
depth. 

[0023] In a .preferred embodiment the acquisition of the 
fluorescence signal is performed in the reflective mode, 
that is to say, the injection and detection points are 
substantially equidistant from the volume comprising the 
fluorophore . 

[0024] Thus in one embodiment of the invention the depth 
of a volume comprising a fluorophore can be determined by 
establishing a calibration curve for which measurements 
are obtained for known depth of the volume and using this 
calibration curve to estimate the depth of a volume which 
is originally unknown. 
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[0025] It will be appreciated that the depth of more than 
one fluorophore can be assessed by detecting a different 
emission wavelength for each of the fluorophore provided 
that their emission spectra do not overlap. 

[0026] In another aspect of the present invention there is 
provided a method wherein an analytical relationship 
between tnax and optical characteristics of the medium is 
used to calculate the depth of the fluorophore. 

[0027] Under certain assumptions such as assuming that the 
optical properties of the medium are the same at the 
excitation and emission wavelength, the fluorescence 
intensity as a function of time can be expressed by the 
Bora approximation: 



_0 

r 



r 



* (impulseRcsponse) 



(equation 1) 
[0029] Where: 

[0030] r cp ia the distance from source (s) to fluorophore 
depth position (p) ; 

[0031] r p a is the distance from fluorophore depth position 
(p) to detector (d) ; 

[0032] ^ is the optical absorption coefficient; 
[0033] D is the optical diffusion coefficient; and 
[0034] v. is the speed of light in the medium 
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[0035] By setting the first derivative of equation 1 as a 
function of time equal to zero, the maximum of the TPSP 
(tw) can be found. Under certain approximations, and by 
assuming that r ap is approximately equal to rp^it is found 
that: 

[0036] t«m = -^=£ (equation 2) 

_ 1 

[0037] where d is the depth of the fluorophore, Dst ~$£' 

It will be appreciated that t ra «, or any time point between 
the earliest detected signal (time »0») up to and including 
tnax (and therefore d) can be more . accurately determined 
numerically using the full equation 1. 

[0038] As can be seen the depth of the fluorophore is 
proportional to t^. In order to determine the depth of a 
volume comprising a fluorophore in a turbid medium, the 
speed of light and the scatter coefficient can be provided 
using known values for the medium. Alternatively, the 
scattering coefficent can be provided using the standard 
time domain approach by direct measurements. 

[0039] The scattering and absorption coefficients can be 
determined at either the excitation and/or emission 
wavelengths of the fluorophore, since in a preferred 
embodiment the method assumes chat these values are the 
same for the excitation and emission wavelengths. However, 
it is also possible to derive t« by using the scattering 
and absorption coefficients as determined at both the 
excitation and emission wavelengths of the fluorophore. In 
this respect, it will be appreciated that a more general 
form of equation l in which the absorption and scatter 
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coefficients, determined at both the emission and 
excitation wavelengths, are -used, is also encompassed in the 
present invention. 

[0040] The calculation is based in part on the assumption 
that the distance between the point of light injection at 
the surface and the f luorophore and between the f luorophore 
and the point of fluorescence emission detection at the 
same surface is substantially the same (figure 1) . 
Determination of the injection and detection points 
relative to the volume of interest for this reflection 
geometry can be achieved, for example by obtaining a 
topographic image of the fluorophore in a region of 
interest. In this respect the 2D (topographic) image may be 
obtained by optical modalities such as TD, CW and frequency 
domain (FD) . 

[0041] Once the injection and detection points have been 
determined, TPSF data is then acquired by injecting pulses 
of light at an excitation wavelength of the fluorophore at 
the injection point and by detecting as a function of time 
the fluorescence emission at an emission wavelength. TPSF 
generation using TD optical method are well known in the 
art. The time at which the TPSF reaches a maximum (t^) is 
then recorded and used in the above relationships (either 
equation 1 or equation 2) to determine the depth of the 
fluorophore. 

[0042] It will be appreciated that a single 
injection/ detect ion points pair is sufficient to determine 
the depth, therefore providing a rapid method for 
determination of depth. However, it is also possible to use 
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multiple injection/detection pairs to improve the accuracy 
of depth determination. 

[0043] For a given depth, the surface CW intensity can be 
related to fluorophore concentration by the optical 
properties of the medium (absorption and scattering 
coefficients) . Thus, in another aspect of the invention, 
estimates of the relative concentration of the fluorophore, 
Cone. Relative* can be obtained by determining its depth, d, 
using the method described above, and normalizing the 
surface CW intensity measurement, CW, as follows (Equation 
3) : 



[0044] 



(equation 3) 



[0045] under certain assumptions, Equation 3 
derived from equation 1. 



can be 



[0046] This method provide a means for rapidly compare the 
relative concentration of a fluorophore at 2 or more 
distinct locations. The determination of the relative 
concentration of- a f luorescently labeled antibody in 
different organs of an animal provides a non- limiting 
example. 



[0047] Relative concent ration! 
estimates of absolute 
measurement is absolutely 
calibration curve can be est 
concentration of the 
abosrption coefficent can be 
for the medium or alternatively 
standard time domain approach 
will also be appreciated that 



can be further processed to 
if the CW surface 
calibrated. That is to say a 
ablished for CW using known 
ore. The scattering and 
provided using known values 
may be provided using the 
by direct measurements. It 
the CW surface intensity may 



concentration 



fluoroph 
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be provided by temporally integrating the time domain 
fluorescent measurement. 

[0048] By further analyzing the CW fluorescence intensity 
data from the above experiments. it was possible to 
calibrate the optical attenuation of the turbid medium. In 
other words. for the inclusion of fixed fluorophore 
concentration, the optical attenuation can be obtained by 
fitting the decrease of CW fluorescence intensity measured 
by the surface detector- as the inclusion depth increases. 
The same optical attenuation can be used to describe the 
decrease of CW fluorescence intensity of fluorophore with 
other concentrations as the inclusion depth increases. 

[0049] Given the above series of experiment . the inclusion 
may be submerged at unknown depth and concentration in the 
turbid medium. As stated above, direct cw fluorescence 
intensity measurements alone cannot decouple the inclusion 
depth from its fluorophore concentration. However TD 
experiments can provide t majt which yields the inclusion 
depth. Once the inclusion depth is known, the calibrated 
optical attenuation of the turbid medium can be used to 
normalise the CW fluorescence intensity measurement to 
yield the relative concentration of the fluorophore 
inclusion. Further, an absolute calibration of the CW 
intensity measurements permits estimates of the absolute 
concentration . 

[0050] In yet a further embodiment of the invention there 
» is provided a method for generating a tomographic image of 
a fluorophore distribution. Topographic images of a region 
of interest can be obtained using well known method such as 
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CW or TD optical imaging. Such topographic images can be 
combined with the depth and/or concentration information 
obtained by the methods described above to generate a 3- 
dimension (tomographic) image of a volume of interest 
comprising one or more f luorophores. 

[0051] It is also obvious to those skilled in the art that 
the above time domain approach can be obtained with 
frequency domain approach via the Fourier transform 
relationship. 

[0052] The embodiment (s) of the invention described above 
is (are) intended to be exemplary only. The scope of the 
invention is therefore intended to be limited solely by the 
scope of the appended claims. 
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I/WE CLAIMs 



A method for determining depth of a volume comprising 
a fluorophore in a turbid medium using time 
domain (TD) optical fluorescence, said method 
comprising: 

i) obtaining Temporal Point Spread Function (TPSF) 
data corresponding to at least one TPSF by 
injecting light at an injection point at an 
excitation wavelength of said fluorophore and 
detecting light at a detection point at an 
emission wavelength of said fluorophore and 
wherein said injection and detection points are 
in a reflection geometry and are substantially 
equidistant from said volume; 

ii) determining a time (t^) substantially 
corresponding to the maximum of said TPSF; 

iii) correlating said tmax with said depth, wherein 
said depth is insensitive to fluorophore 
concentration . 

2. The method as claimed in claim 1 wherein said step of 
correlating comprises: 

a) establishing a calibration curve of t m3 x vs depth 

for a plurality of depths; 

b) using said calibration curve and said tmax 
determined in step ii) to obtain said depth of 
said volume. 

3. The method as claimed in claim 1 wherein said step of 
correlating comprises: 
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a) providing at least one optical property of said 

medium, speed of light in said medium and a 
lifetime of said f luorophore; 

b) determining a calculated depth using said tmax and 

said at least one optical property, said speed of 
light and said lifetime. 

4. The method as claimed in claim 3 wherein said at 
least one optical property is scatter coefficient. 

5. The method as claimed in claim 1 further comprising a 
step of estimating a position of said fluorophore 
prior to said step of obtaining said TPSF. 

6. The method as claimed in claim 5 wherein said step of 
estimating is performed by obtaining a topographic 
image of a region of interest containing said 
fluorophore . 

7. The method according to claim 1 wherein said 
fluorophore is an intrinsic fluorophore. 

8. The method as claimed in claim 3 wherein said at 
least one optical property is obtained using time 
domain optical measurements of said medium. 

9. The method as claimed in claim 3 wherein said at 
least one optical property is estimated. 

10 . The method as claimed ■ in claim 3 wherein said at 
least one optical property is obtained by matching 
said medium with predetermined optical properties of 
media in a database. 
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11. 



12. 



13 • 



14. 



15. 



The method as claimed in claim 3, wherein Baid at 
least one optical property and said speed of light 
are substantially the same at both the excitation or 
emission wavelength of the fluorophore and are 
determined at either said emission or said excitation 
wavelength. 

The method as claimed in claim 3, wherein said at 
least one optical property and said speed of light 
are determined at the fluorophore excitation and 
emission wavelength. 

The method as claimed in claim 3, wherein said at 
least one optical property and said speed of light 
are determined taking into account optical properties 
of the turbid medium at said fluorophore location. 

A method for estimating concentration of fluorophore 
in a volume in a turbid medium using optical 
fluorescence, said method comprising: 

i) obtaining depth of said volume; 

ii) providing optical properties for said medium; 

iii) obtaining a CW intensity surface reflection 
measurement of said fluorophore; and 

iv) normalizing said CW intensity measurements using 
said optical properties and said fluorophore 
depth to obtain a relative fluorophore 
concentration . 

The method as claimed in claim 14 wherein said depth 
is obtained by the method of any one of claim 1-13. 
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The method as claimed in claim 14 wherein said 
optical properties are obtained using time domain 
optical measurements of said medium. 

The method as claimed in claim 14 wherein said 
opcical propercies are estimated. 

The method as claimed in claim 14 wherein said 
optical properties are obtained by matching said 
medium with predetermined optical properties of media 
in a database. 

19. The method as claimed in claim 14, wherein said 
optical properties are substantially the same at both 
the excitation or emission wavelength of the 
fluorophore and are determined at either said 
emission or said excitation wavelength. 



17 



18. 



20. 



The method as claimed in claim 14, wherein said 
optical properties are determined at the fluorophore 
excitation and emission wavelength. 

i 

21. The method as claimed in claim 14, wherein said 
optical .properties are determined taking into account 
optical propercies of the turbid medium at said 
fluorophore location. 

22 . A method for determining concentration of fluorophore 
in a volume in a turbid medium using optical 
fluorescence, Baid method comprising: 

i) obtaining depth of said volume; 

ii) providing optical properties for said medium; 
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iii) obtaining a CW intensity surface reflection 
measurement of said fluorophore; and 

iv) calibrating said CW intensity surface reflection 
measurement so as to obtain an absolute 
concentration of said fluorophore. 

23. The method as claimed in claim 22 wherein said depth 
is obtained by the method of any one of claim 1-13 . 

24. The method as claimed in claim 22 wherein said 
optical properties are obtained using time domain 
optical measurements of said medium. 

25. The method as claimed in claim 22 wherein said 
optical properties are" estimated. 

26. The method as claimed in claim 22 wherein said 
optical properties are obtained by matching said 
medium with predetermined optical properties of media 
in a database. 

27. The method as claimed in claim 22, wherein said 
optical properties are substantially the same at both 
the excitation or emission wavelength of the 
fluorophore and are determined at either said 
emission or said excitation wavelength. . 

28. The method as claimed in claim 22, wherein said 
optical properties are determined at the fluorophore 
excitation and emission wavelength. 
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The method as claimed in claim 22, wherein said 
optical properties are determined taking into account 
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optical properties of the turbid medium at said 
fluorophore location. 

30. A method for generating a tomographic image of a 
fluorophore distribution in a turbid medium said 
method comprising: 

i) obtaining a topographic image of said fluorophore 

distribution; 

ii) determining depth of a plurality of volumes of 

interest comprising said fluorophore using the 
method as claimed in any one of claim 1-13; 

iii) combining said depth information and said 
topographic image to generate a tomographic image 
of said distribution. 



The method as claimed in claim 30 wherein said 
tomographic image is further processed with the 
method as claimed in any one of claim 14-29 to 
generate a tomographic relative flurorophore 
concentration or absolute fluorpohore concentration 
image . 

32 . A method for determining a relative or absolute 
concentration of a fluorophore in a turbid medium 
said method comprising: 

i) establishing a calibration curve relating a 

continuous wave intensity measurement in said 
turbid medium and concentration and depth of said 
fluorophore; 

ii) determining a depth of said fluorophore using the 

method as claimed in any one of claim 1-13; and 
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iii) determining said concentration using said 
calibration curve. 

3. A method for determining depth of a volume comprising 
a fluorophore in a turbid medium using time 
domain (TD) optical fluorescence, said method 
comprising: 

i) providing an equation relating fluorescence and 
time; 

ii) obtaining Temporal Point Spread Function (TPSF) 
data corresponding to at least one TPSF by 
injecting light at an injection point at an 
excitation wavelength of said fluorophore and 
detecting light at a detection point at an 
emission wavelength of said fluorophore and 
wherein said injection and detection points are 
in a reflection geometry and are substantially 
equidistant from said volume; 

iii) numerically solving said equation using said 
obtained TPSF data to determine a predetermined 
time t comprised between 0 and t ma x;and 

iv) correlating said time t with said depth, wherein 
said depth is insensitive to fluorophore 
concentration. 

34. The method as claimed in any one of claims 1-13, IS, 
16, 23, 24 and 30-33 wherein time domain informati6n 
is obtained by acquiring Frequency Domain (FD) data 
and applying Fourier transform to said data. 
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